Mon. Not. R. Astron. Soc. 000, [T]-?? (2008) Printed 8 June 2009 



(MN MfejX style file v2.2) 



Pixel- lensing as a way to detect extrasolar planets in M31 

G. Ingrosso/* S. Calchi Novati,^ F. De Paolis,^ 
b^'Ph. Jetzer,^ A. A. Nucita^ and A.F. Zakharov^'^ 

C " 3 ' ^ Dipartimento di Fisica, Universita del Salento and INFN Sezione di Lecce, CP 193, 1-73100 Lecce, Italy 
C 2 ^ ^ Dipartimento di Fisica, Universita di Salerno, 1-84081 Baronissi (SA) and INFN Sezione di Napoli, Italy 
^Nj , ^ Institute for Theoretical Physics, University of Ziirich, Winterthurerstrasse 190, CH-8057 Ziirich, Switzerland 
* XMM-Newton Science Operations Centre, ESAC, ESA, PO Box 50727, 28080 Madrid, Spain 
^ Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, 117259 Moscovi, Russia 
^ Bogoliubov Laboratory of Theoretical Physics, loint Institute for Nuclear Research, 141980 Dubna, Russia 



Accepted xxx; Received xxx; in original form xxx 



ABSTRACT 



We study the possibility to detect extrasolar planets in M31 through pixel- lensing 
observations. Using a Monte Carlo approach, we select the physical parameters of the 
binary lens system, a star hosting a planet, and we calculate the pixel-lensing light 
curve taking into account the finite source effects. Indeed, their inclusion is crucial 
since the sources in M3f microlensing events are mainly giant stars. Light curves with 
detectable planetary features are selected by looking for significant deviations from the 
corresponding Paczyhski shapes. We find that the time scale of planetary deviations 
in light curves increase (up to 3-4 days) as the source size increases. This means that 
only few exposures per day, depending also on the required accuracy, may be sufhcicnt 
to reveal in the light curve a planetary companion. Although the mean planet mass 
for the selected events is about 2 Mjupiter, even small mass planets [Mp < 20 Mq) 
can cause significant deviations, at least in the observations with large telescopes. 
However, even in the former case, the probability to find detectable planetary features 
in pixel-lensing light curves is at most a few percent of the detectable events, and 
therefore many events have to be collected in order to detect an extrasolar planet in 
M31. Our analysis also supports the claim that the anomaly found in the candidate 
event PA-99-N2 towards M31 can be explained by a companion object orbiting the 
lens star. 
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1 INTRODUCTION 



In the last years, it has become clear that gravitational mi- 
crolensing, initially developed to search for MACHOs in our 
Galactic halo and near the Galactic disk l|Paczvriskilll986l: 
Alcock et ail Il993l: iPaczvriskil Il996l : iRoulet fc MollerachI 



19971 . l2002l : IZakharov fc SazhinHl99i) can be used to infer 



the presence of extr asolar planets orbiting around lens stars 
(see the revie w bv lPerrvmanll2000l : IPerrvman et Z ll2005l: 
lBennett|[2009l '). 

As shown bv lMao fc Paczvriskil (|l99ll ) the planet pres- 
ence effect on the light curve in a microlensing event to- 
wards the Galactic bulge is generally a short duration per- 
turbation to the standard microlensing curve. These de- 
viations last from a few hours to some days (depending 
on the planet mass) and can occur relatively frequently. 
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even for rather small mass planets. Indeed, the microlensing 
technique is sensitive to planets in a rather large range of 
masses, from Jupiter-l ike planets down to Earth-like ones 
(|Bennett fc Rhielll996D . 

iGould fc Loebl (|l992) pointed out that there is a signifi- 
cant probability to detect planets around stars in the Galac- 
tic disk that act as microlenses by magnifying the light of 
observed stars in the Galactic bulge. Until now, the detec- 
tion of eight extrasolar planets has been reported by using 
the microlensing technique (IBond et al . 2004; Udals ki et al. 
|2005j : [Beau Ueu et"al] l2006l : iGould et al.. .2006. : .Gaudi et al. 
2008;:; ^Bennett 20oi). We remind that the masses of three 
of them (~ 3, 5 and 13 M^) are at the lower bound of the 
detected planetary mass range. Indeed, more than 300 extra- 
solar planets discovered until now by radial velocity, transit 
and direct imaging m ethods are biase d towards large mass 
(Jupiter- like) planets l|lda fc Linll2004l ). However, radial ve- 
locity searches by ground based experiments have now pro- 
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vided extrasolar planets with Mmin = 2 (|Mavor et alj 
I2OO9I ). whereas space based observations are expected to de- 
tect many Earth-mass planets (with Kepler satellite Q) and 
many Earth-size planets (with COROT spececraft . 

A further advantage of the microlensing is that it 
works better for large distance of the source star, since 
the optical depth increases by increasin g the d istance, as 
one can already see from the lEinsteinl (1936) approach. 
This gives the opportunity to detect planetary systems 
at distances much larger with respect to those accessi- 
ble by the other te cniques and even i n other galaxies 
such as M31 (see, e.g. JCovone et al1l2000l : iBaltz fc Gondold 
I2OOII ). In this case, however, the source stars are not re- 
solved by ground based telescopes - a situation referred to 
as "p ixel-lensing" l|Crottsl [1993 : iBaillon etHI Il993l : iGouldl 
1 19961 ) - and only bright sources (i.e. giant stars with large 
radii), sufficiently magnified, can give rise to detectable mi- 
crolensing events l|Ansari et al.|[l997l ). This implies that fi- 
nite size effects, leading to smaller planetary deviations in 
pixel-lensing light-curves with respect to microlensing to- 
wards th e galactic bulge, cannot be neglected (see, e.g., 
iRiffeser. Seitz fc Bende r 2008). Usually, highly magnified 
events arise when the source and lens stars align very 
closely. In this case there is the largest chance of observ- 
ing the perturbations in the light curves induced by plan- 
ets l|Griest fc Safizadeh|[l998h . This is particularly true for 
large mass planets, for which the planetary signals are not 
strongly suppressed by finite size effects, whereas for low 
mass planets, the planetary sign als may remain detectable 
during other phases of the event ()Bennettll2009t ). 

Until now, only about a dozen microlensing events 
have been observed t oward s M31 by the POINT- AGAPE 
Calchi Novati et al, l2005f ) and MEGA collaborations 
de Jong et al Only in one case a deviation from 



the standard Paczyhski shape has been observed and at- 
tributed to a se condary component orbiting the lens star 
ijAn et all l2004l ). However, new obser vational campaigns 
towards M31 have been und ertaken (jKerins et al.l I2OO6I : 
ICalchi Novati et al.|[2007l , |2009| ) and hopefully a few planets 
might be detected in the future, providing a better statistics 
on the masses and orbital radii of extrasolar planets. It is in 
fact expected, and supported by observations and numerical 
simulations, that almost a ny star has at least a planet or- 
biting around it (see, e.g.. iLinew eaver fc Grethcr 2003). In 
other words, as also suggested by Baltz_fc Gondolo (200l|), 
the rate of single lens events towards M31 may suffer of a 
strong contamination of binary lensing events, most of which 
are expected to be due to extrasolar planets. 

Therefore, it is important to address the question of how 
to extract information about planetary lensing events, oc- 
curring in M31, from the observed microlensing light curves. 
Since planetary perturbations last from hours to a few days, 
a monitoring program with suitable sampling must be re- 
alised, in order to avoid missing these perturbations. The 
feasibility of such rese a rch p rogra m has bee n alrea dy ex- 
plored bv IChung et~al] lj2006l ) and iKim et aD l|2007l ). They 
have considered the possibility to detect planets in M31 
bulge by using the observations taken from the Angstrom 
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collaboration (jKerins et al.|[200^ with a global network of 
2 m class telescopes and a monitoring frequency of about 
five observations per day. The analysis for planet detection, 
however, has been performed by using a fixed configuration 
of the underlying Paczyhski light curve. 

In the prese nt work, instead, by using a Monte Carlo 
(MC) approach l|De Paolis et al.ll2005l : Ungrosso et al.ll2006l . 
|2007f) we explore the possibility of detecting extrasolar plan- 
ets in pixel-lensing observations towards M31, by consider- 
ing the multi-dimensional space of parameters for both lens- 
ing and planetary systems. Taking into account the finite 
source effects and the limb darkening and using the residual 
method we can select the simulated light curves that show 
significant deviations with respect to a Paczyhski like light 
curve, modified by finite source effects. The advantage of the 
Monte Carlo approach is that of allowing us a complete char- 
acterization of the sample of microlensing events for which 
the planetary deviations are more likely to be detected. 

The paper is structured as follows. In Section 2, we 
give the basics of binary-lensing events. In Section 3 we dis- 
cuss our MC simulations for planetary detection in M31. In 
section 4 we present our main results and in Section 5 we 
address the conclusions. 



2 BINARY-LENSING EVENTS 
2.1 Generalities 

If a source star is gravitationally lensed by a binary lens, 
the equation of lens mapping from the lens plane t o the 
source plane can be exp ressed in complex notation l|Wittl 



source plane can be exp 
ll990l : IWitt fc Maolll995r ) 



E 



TUj/M 
Z - ZL,j 



(1) 



where — ( + irj and z = x + iy are the source and the 
image positions, 2 is the complex conjugate of z, mi, m2, 
and zl,2 are the the masses and the positions of the 
two lenses, respectively. Here and in the following, all the 
lengths (angular separations) are normalized to the radius 
-Re (angle ^e) of the Einstein ring which are related to the 
physical parameter of the lens by 



-Re = 



/ 4GM \ Di^{Ds-Dl) 
V c2 J -Ds 



1/2 



and 



(2) 



where M — mi + m2 is the total mass of the binary system, 
Dl and -Ds are the distances to the lens and to the source, 
respectively. Under the condition mi > m2, we define the 
mass ratio parameter q = m2/mi. In addition, we assume 
that the two masses of the binary system are located on the 
real axis, with the centre of mass in the origin. Let us denote 
with d the angular separation between the two objects in 
units of 6e- 

To determine the image position and magnification, one 
has to take the complex conjugate of equation ((TJ and sub- 
stitute the expression for z back in it, obtaining a fifth-order 
polynomial in 2, i. e. p{z) = X]^_i CiZ^ = (with coefficients 
Ci depending on M, d, and g), whose solutions give the im- 
age positions. Due to lensing, the source star image splits 
into several fragments up to a total number Nj. Since the 
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lensing process conserves the source brightness and thus the 
magnification of each image, the tot al magnification co rre- 
sponds to the sum over all images (jWitt fc Madlf995l ). i. 



N, 



( Pjzi) 
\AeiJ ) ' 



(3) 



where the determinant of the Jacobian is 



det J : 



dz dz 



(4) 



A planetary lens system is characterized by the condi- 
tion that the planet mass Mp = m2 is much smaller with re- 
spect to the host star mass Ml = mi. fn this case, the planet 
only induces a perturbation on the underlying Paczyiiski 
curve of the primary lens. Planet perturbations occur when 
the source star crosses and/or passes near caustics, which 
are the set of source positions on the 77) plane at which 
the magnification is infinite (i.e. those corresponding to 
det J = 0) in the idealized case of a point source. Clearly, for 
realistic sources of finite size the magnitude gets still quite 
large, but finite l|Witt fc Mao f994). Caustics form a single 
or multiple sets of closed a nd concave curves (fold caustics) 
which meet in cusp points llSchneider. Ehlers fc Falcd 1 19921 : 
[Schneider fc Weisslll992l : IZakharov fc Sazhinlll995f rThe lo- 
cation of the planet perturbations depends on the position 
of the caustics and the source trajectory. 

There have been several attempts to determine caus- 
tic positions and shapes by using analytic methods and 
treating the planet induced deviations as a per turbation 
jGaudi fc Gouldlll997l : [bominik|[l999l : lBozzalll999l ). For the 
planetary case, there exists two sets of caustics: "central" 
and "planetary". The single, central caustic is located on 
the star-to-planet axis, close to the host star. For a wide 
range of parameters the caustic has a diamond shape and 
ca n be described by paramet ric equations (as it was shown 



bv lZakharov fc Sazhinlll997a 



if the lChang fc Refsdall 



1984, 



, central astroid caustics arise 
model is used). Planetary 
caustics are located away from the host star, at distance 
~ {cP — l)/d from the primary lens position. There is one 
planetary caustic (with a diamond shape) on the star-to- 
planet axis, on the planet side, when d > 1 and two sets 
of caustics, off the axis, (with triangular shape) on the 
star side when d < 1. The dimensions of both central and 
planetary caustics increas e by increasi n g the mass ratio q 
JZakharov fc SazhinI [l99i : lBozzalll999l : IChung et allbooa 
iHan fc Gaudill2008l ). Moreover, for a given q value, the caus- 
tic sizes are maximized when the planet is inside the so 
called "lensing zone" , which is defined (with some arbitrari- 
ness) as the range of star-to-planet separation 0.6 <d< 1.6 

l|Gould fc Loebl Il992l : [ Griest fc Safizadeh|[l998l ). The time 
duration scale of the perturbations induced by a planet 
and the probability of their detection are proportional to 
the caustic size, at least when this region is large enough 
so that the planet ary signals are not suppressed by the fi- 
nite source effects llMao fc Paczviiskilll99ll : iBolatto fc Falcd 
ll994l : lGould fc Loeblll99^ 



2.2 Finite source approximation 

Since in pixel-lensing towards M31 the bulk of the source 
stars are red giants (see Section 3), one has to take into 
account the source finiteness. This leads to smaller plane- 
tary deviations in pixel-lensing light curves with respect to 
microlensing towards the galactic bulge, for which the point- 
like source approximation is acceptable. For finite source ef- 
fects with limb darkening t he magnification has to be numer- 
ically evaluated (see, e.g ., Schneider, Ehlers fc Falcd [199I : 
iBoed anov fc ChereDashchukiil995a : iDominikii2005r and ref- 
erences therein) 



CdeJ^'A^{e,~p;t)i{p)pdp 

2n J^'l{p)pdp 



(5) 



where p = 0s/9e is the normalized angular size of the source 
{9s ~ Rs/Ds, Rs being the source radius), and I{p) is the 
intensity profile of the source i ncluding limb darkening, for 
which we use the lClared l|2000l ') approximation 

I{p) = l-ai(l-/ii/'')-a2(l-M)-a3(l-M'''')-a4(l-M') ,(6) 

with p = p/p and the coefficients in the 7?- band are ai — 
0.8282, a2 = -0.9866, as = 1.6801 and 04 = -0.6604 (for 
red giant stars). 

Moreover, since during caustic crossing the magnifica- 
tion could have strong changes and (at least for small mass 
planets and/or realistic source sizes) typical time scale for 
crossing could be comparable with the exposure time tcxp 
(needed to have a reasonable signal-to-noise level) we take 
the average magnification of equation ([Sjl in the interval 

(^ icxp/2, t tcxp/2). 

Finite size source efi^ects can be relevant for two reasons. 
First, the relationship between the dimensionless radius p 
and the impact parameter uo determines if the finite size 
effects are important or not for the main microlensing light 
curve. This occurs in the events with p/uo > 1 or p/uo < 1, 
respectively. Second, finite size effects may be important for 
the planetary deviations even if they are not relevant for 
microlensing without planets. Indeed, p is to be compared 
not only with uo, but also with the caustic size A. In partic- 
ular, whenever p/uo > 1, it results that p is typically much 
larger than A (at least for small enough mass planets). In 
this case, smoothed planetary deviations are produced in 
the light curves, since the planetary magnification has to be 
averaged on the source area. In a similar way, depending on 
the lens system geometry and proper motion, whenever the 
ratio p/uo < 1, stronger and temporally localized planetary 
deviations are produced in the light curves, since the caustic 
region results to be a non negligible fraction of the source 
area. Within the following analysis for the detection of plan- 
etary deviations we are going to identify two classes (I and 
II) of events, depending on the ratio p/uo > 1 and p/uo < 1, 
respectively Q. 



http:/ /webviz.u-strasbg.fr/viz-bin/VizieR- 
source=J/A-f A/363/1081. 

^ We mention the cla ssification of the planetary perturbations 
by ICovone et al. I 1 I2OOOI) . which distinguishes two main types of 
anomalies in the light curves, namely the events aflpected by the 
central caustic (type I) , and the ones aff'ected by one of the plan- 
etary caustics (type II). In our analysis we do not attempt to 
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3 MONTE CARLO SIMULATION 
3.1 Light curve generation 

In the present analysis we assume that the lens is a binary 
system constituted by a star and a planet companion 0. Our 
aim is to evaluate the probability to detect the presence of 
planets in M31 through Earth-based pixel- lensing observa- 
tions with telescopes of different diameters. These telescopes 
could be initiated to observe towards a microlensing event 
candidate, so making a high cadence observations of an on- 
going microlensing event. As reference values, we adopt a 
CCD pixel field of view of 0.2 arcsec, a typical seeing value 
of 1 arcsec and an average background luminosity at tele- 
scope site of ~ 21 mag arcsec"^ in ii-band. To have a good 
S/N ratio we consider in the MC analysis exposure times 
tcxp of 30 minutes. Moreover, we assume a regular sampling 
neglecting any loss of coverage due to bad weather condi- 
tions. 

In order to take into account the spatial variation 
of the background level we select four directions (named 
A, B, C, D) at increasing distances from the M31 centre. 
Assuming a coordinate system with origin in the M31 centre 
and axes along the north-south and east- west directions, the 
coordinates of the selected directions are the following: A (- 
6,0) arcmin, B (-9,0) arcmin, C (-12,0) arcmin, D (-21,-6) 
arcmin. In the direction A the microlensing is primarily due 
to self-lensing events by stars in the M31 bulge and disk, 
whereas towards the external directions the contribution to 
microlensing due to lenses belonging to the M31 halo be- 
comes larger. Our investigation of the D direction is moti- 
vated by the det ection of the a nomaly in the pixel-lensing 
event PA-99-N2 jAn et al.ll20ol ). 

As for the generation of the trial microlensing 
light curves we cl osely follow the approach outlined by 
iKerins et al.l ll200l|). The adopte d M31 astrophysical model 
was described by Ingrosso et al.l l|2006l ). Once the event lo- 
cation has been selected, for any lens and source population 
lying along the line of sight, we use a MC approach to select 
the physical parameters of the systems: source magnitude, 
primary lens mass, source and lens distances, effective trans- 
verse velocity of source and lens, impact parameter of the 
lens. 

The luminosity of star sources, mainly red giants in 
the interval of absolute magnitude (—4, 2.4), and the corre- 
sponding radii are drawn from a sample of stars generated 
by a synthetic color -magnitude diagram comp utation algo- 
rithm|j described bv lAparicio fc Gallarj (|2004|) based on the 
stellar evolution library ( Bertelli et al.lll994h a nd the bolo- 
metric correction database ( Girardi et al.ll2002f ). 

As next, we have to select the mass Mp and the (pro- 
jected) orbital distance dp of the extrasolar planet. Most 
of the hundreds of extrasolar planets discovered up to now 
(see the web site http://exoplanet.euJ have typically very 



characterize the planet deviations as due to the intersection of 
central and/or planetary caustics, but we look at the shape of 
the induced planetary features on the light curves. A classifica- 
tion of the events based on the caustic crossing will be considered 
in a forthcoming paper. 

^ Ba sed on the recent detections of multiple planets jGaudi et al.l 
l2008t) . one can expect that this assumption is rather conservative. 
^ http:/ /iac-star. iac.es/iac-star. 



large mass es and orbit at small distances around their par- 
ent stars l|Udrv fc Santod 120071). This appea rs to be a re- 
sult of observational biases Jida fc Linll2004h since most of 
the planets have been detected by radial velocity and tran- 
sit techniques that are most sensitive to massive and close 
planets. Direct imaging and microlensing techniques con- 
tribute only a minor fraction of the detected events. Indeed, 
available theoretical and numerical analysis show that most 
extrasolar planets are expected to have relatively smaller 
masses. Furthermore, the (projected) orbital distance from 
their host ing stars is expected in the range ~ 0.04— 100 AU 
(see, e.g.. lTabachnik fc Tremainell2002l : llda fc Linll2004h . In 
the present paper we assume that the distribution of Mp 
and orbital period P, for M p < lOMjupiter, is given by th e 
simple analytical expression ([Tabachnik fc Tremainell2002l ) 



dn{Mp,P) = C Mp" P 



( 



dMp 



V Mp 



(7) 



with a — 0.11 and /3 — —0.27. This relation is obtained by 
investigating the distribution of masses and orbital periods 
of 72 extrasolar planets, taking into account the selection 
effects caused by the limited velocity precision and duration 
of existing surveys. We note that in the analysis leading 
to the above distribution, it was assumed that the stars in 
the survey are of solar type, and therefore any dependence 
(as implied by recent extrasolar planet observations) of the 
planet mass on the parent star mass and metallicity has been 
neglected. Taking that into account, one would certainly re- 
place equation (O with a different one, and therefore the 
results presented in Section 4 for the detectable planet rate 
would change. For example, a steeper planet mass distri- 
bution ( a s tha t found for all Doppler-detected planets by 
IJohnsonI l|2009l ) with a = 0.4) implies a smaller (about a 
quarter) overall planet detection rate, as a consequence of 
the decrease of the mean planet mass. More importantly, 
a dependence of the planet mass distribution on the par- 
ent star mass would introduce a dependence of the planet 
detection rate on the lens population (bulge or disk stars) 
that could be recognized, provided a sufficient event statis- 
tics towards different lines of sight would be available. In 
equation ((7)l , the upper limit of the planetary mass is set at 
Alp = 10 Mjupiter- This roughly corresponds to the usually 
assumed lower mass limit for brown dwarfs. Indeed, in the 
range 10 — 20 Mjupitor the two populations overlap. More- 
over, in the simulation we select a lower planetary mass limit 
of 0.1 M^. Once the masses of the binary components and 
the planet period have been selected, the binary separation 
dp is obtained by assuming a circular motion of the planet. 

As a parameter in our MC analysis we introduce the 
number Nim of images per day. We take Nim in the range 
2-12 day~^, the latter value corresponding to a sampling 
time of two hours. For all selected values of A^im, the corre- 
sponding binary light curve at any time is given by 



5p(t) = /bi + /o((Ap(t))-l) 



(8) 



where /bi is the background signal from the galaxy and the 
sky, /o is the unamplified source star flux and {Ap{t)) the 
time varying magnification, that takes into account both the 
source finiteness and the motion of the lens-source-observer 
system duri ng the exposure tim e texp. To mimic superpixel 
photometry l|Ansari et al.|[l997h used in a real observational 
campaign we evaluate the star and the background flux 
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within a n-pixel square "superpixel" , whose size n is de- 
termined to cover most of the average seeing disk. We recall 
that we consider the pixel-lensing regime where the noise 
is do minated by the (lin e of sight dependent) background 
noise l|Kerins et alJlioOll ). Accordingly, we add to Sp(t) a 
Gaussian noise. 



3.2 Microlensing event selection 

As a first step, within the MC simulation, we have to test 
whether the flux variation due to the microlensing event 
is significant with respect to the background noise a{x,y), 
where {x, y) identifies the line-of-sight. To asses the detec- 
tion of a flux variation we evaluate its statistical significance 
testing whenever and to what extent at least three consec- 
utive points exceed the baseline level by 3ff, f ollowing the 
analysis described bv lCalchi Novati et al.l()2002l ). We remark 
that the condition on the variation significance is the only 
one used at this stage. In the following we refer to events that 
show a significant fiux variations as to "detectable" events. 



3.3 Planet detection 

The expected signature of an extrasolar planet orbiting the 
lens star is the presence of perturbations with respect to the 
corresponding smooth Paczyhski light curve. Therefore, we 
look for a selection criterion based on the analysis of the sig- 
nificance of such deviations. To this purpose, given the wide 
range of the microlensing parameters and the corresponding 
planetary deviations, we consider two indicators for which 
we select (by the direct survey of many light curves) thresh- 
old values. They are the mean deviation (in units of a) of the 
planetary light curve from that of a single lens event, and 
the maximum value of the time dependent relative planetary 
magnification (in units of the expected Paczynski value) . 

At first, we fit the light curve in equation ((8)l with a 
Paczynski-like law modified to take into account finite source 
effects and determine the best fit parameters. The latter 
are the baseline flux /^i, the maximum magnification time 
to, the unamplified star flux /o, the Einstein time the 
dimensionless impact parameter Mq and the dimensionless, 
projected star radius p". Accordingly, the time dependent 
flux S''\t) due to a single lens event is given by 

S°(t) = /^i + /o°({A"(i))-l) , (9) 
where the magniflcation l|Einsteinll 19361 : IPaczvnskill 19861 ) 



A\t) = 



(10) 



is given in terms of the time varying (normalized) lens an- 
gular distance to the source u^{t) 



(11) 



and (v4°(t)) is the analogous of equation ((Sjl, in this case 
evaluated by using the analytical approximation given by 
IWitt fc Maol (f 1994 ). Then, we can evaluate the time depen- 
dent variable 
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Figure 1. Normalized (to unity) distributions of 4^/2 (the du- 
ration of a microlensing event without planet). Top panel: de- 
tectable events. Bottom panel: selected events (xr > 4, Afgooj > 3 
and (e)niax > 0.1) for I class {p/ug > 1, solid line) and II class 
{p/uQ < 1, dashed line) events. In Figs. 1 — 11 we take D = 8 m 
telescope parameters and Mm = 12 day~^. 



and the residual to the single lens fltQ 

x?(t) = (i-x'W)' , 



(13) 



(12) 



where Sp{t) is the light curve including the planet pertur- 
bations, S°(t) the Paczynski flt as abo ve and a{t) is eval- 
uated according to iKerins et all l|200ll ). Therefore, we can 
consider large values of Xr(i) as a significant indicator of 
the presence of detectable planetary deviations in the light 
curves. Actually, we use the sum of the residuals along the 
whole light curve, namely x? = '^iXr{ii)/Ntot, as a first 
quantitative measure of the statistical significance of the 
planetary signals in the ongoing microlensing event. Here 
ti = to + [{t{ — to)/Ntot] i, where to (tf) is the initial (fi- 
nal) instant and A^tot the total number of points. By the 
direct survey of many light curves we select a threshold 
value Xr th = 4. We further require a minimal number 
of points Ngood, even not consecutive, which deviate signif- 
icantly (over 3a) from the Paczyhski best fit. We adopt the 
criterion (i) Xv > th = 4: and A''good > A''good th = 3. 
In other words, the latter condition means that if we have 
significant deviations at only one or two points, we cannot 
conclude that they are caused by a planet orbiting the pri- 
mary lens. 

The light curves fulfilling the above condition (i) may 
show only an overall distortion with respect to the under- 
lying Paczyhsky shape. This is characteristics, in particu- 



^ The analysis of residuals is a well known technique widely ap- 
plied to search for deviations with respect to a null hypothesis. 
* We find that the residuals Xr follow a Gaussian distribution 
(with mean value ~ 1.4 and standard deviation ~ 0.3) in the case 
of light curves generated as single lens events and subsequently 
fitted with a Paczyhski law. 
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Figure 2. Normalized (to unity) distributions of Rmax (the 
magnitude corresponding to the flux variation at the maximal 
Paczynski magnification). Top panel: detectable events. Bottom 
panel: selected events for I class {p/ug > 1, solid line) and II class 
{p/uo < 1, dashed line) events. Here we take D = 8 m telescope 
parameters. 



Figure 3. Event of the I class #1 (see Table [TJ. The upper 
panel shows the simulated light curve (black dots) and the cor- 
responding best fit model (continuous line), that is a Paczynski 
light curve modified for finite source effects. The bottom panel 
gives the difference. 



lar, for events with large source radii and small planetary 
masses. Our second criterion is therefore meant to look for 
and quantify the single more significant planetary pertur- 
bations. To this purpose we consider the time dependent, 
average (with respect to the source area E) relative plane- 
tary magnification 

^^^^^^ ^ r f^d'S[\Ap{5:,t)-A°{x,t)\/A°{x,t)] \ ^^^^ 

V /e '^^^ J 

This quantity is sensibly different from zero only when, de- 
pending on the source and lens parameters and relative mo- 
tion, there is (at a given time) a substantial overlapping 
between the source area and the caustic (central and/or 
planetary) region. So, to select light curves with detectable 
planetary features, besides condition (i), we further require 
that (a) there exist at least one point on the light curve 
with (e)max larger than (e)th ~ 0.1. By using both condi- 
tions the number of selected events get reduced of about 
50% with respect to the events selected by using only the 
condition (i). The condition (i) is particularly efficient to se- 
lect light curves with a large number of points deviating from 
the Paczynski fit, the condition (ii) ensures the presence on 
the light curve of at least one clear planetary feature. Note 
that in this analysis we do not attempt to further charac- 
terize the planet deviations as due to intersection of central 
and/or planetary caustics. 



4 RESULTS 

In the following analysis we consider four different tele- 
scope diameters D = 1.5, 2.5, 4 and 8 m (corresponding 
to zero-point in the _R-band of 23.1 ,24.3, 25.3 and 26.8 
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Figure 4. The same as in Fig. \3\ for the I class event #2 (see 
Table [TJ. 



mag, respectively), toxp = 30 min for all cases and we take 
Nim = 12 day~^, corresponding to a regular sampling time 
of two hours. The effect of taking larger telescopes is that of 
increasing the number of faint, detectable events. Moreover, 
we consider only self-lensing events towards the four consid- 
ered lines of sight (see Section 3.1), leaving out the eventual 
MACHO component in the galactic halos and assume the 
presence of a planet orbiting around each star in the M31 
bulge and disk. 
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Table 1. Parameters of events shown in Figs. [3] -El We also give some microlensing parameters and in the last three columns the sum 
of residuals Xr along the whole light curve, the sum of residuals Xr max and the maximum value of the relative planetary magnification 
{e)max during the time interval corresponding to the strongest planetary feature. 



P/uq Uq dp/R^ Alp 6 /?E -Rmax tl/2 Xr Xr max (e)max 

(Mjupitcr) (dcg) (AU) (day) (mag) (day) 



#1 


2.89 


9.47 X 10-3 


0.90 


4.74 
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2.2 


16.1 
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0.5 
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0.64 


#2 


1.18 


2.63 X 10-2 


0.68 


0.82 
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21.1 


4.0 
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0.12 
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0.57 
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Figure 5. The same as in Fig. [3] for the II class event #3 (see 
Table [TJ . We note that the fit follows the simulated data except 
for a small time interval (5 < t < 10 day). 



An advantage of the Monte Carlo approach to the bi- 
nary microlensing analysis is that we can characterize the 
events with planetary detections. We remind that these 
events have been selected, from the whole sample of de- 
tectable events, by requiring Xi > 4, Ngood > 3 and 
(e)max > 0.1. 

In the Figs. [T] and [2] (for D = 8 m) we give the 
distributions of ti/2 and iZmax for detectable (top panels) 
and selected even ts (bottom panels) 0. As usual (see, e.g., 
iKerins et alll200lh . ti/2 

is the full-width half-maximum mi- 
crolensing event duration and i?max the magnitude in the 
i?-band corresponding to the flux variation at the maximal 
Paczyhski magnification. Comparing the corresponding dis- 
tributions, we see that events with short time duration and 
large flux variation (therefore with smaller impact parame- 
ter) have a larger probability to show planetary deviations. 
This result is due to the fact that the crossing of the central 
caustic (close to the primary lens star) by the source trajec- 



^ We notice that the distributions iJmax for detectable 

and selected events weakly depend on telescope diameter D. 



Figure 6. The same as in Fig. [3] for the II class event #4 (see 
Table [TJ. 

tory is more probable in events with source and lens closely 
aligned. 

As next, for the selected events (bottom panels in the 
Figs. [T] and [2)) we discriminate two classes of events (indi- 
cated by I and II), according to the ratio p/uo > 1 (solid 
lines), or, p/uo < 1 (dashed lines). The ratio p/uo character- 
izes the relative size of the source with respect to the event 
geometry, since p is the dimensionless source radius in the 
lens plane and uo is the dimensionless lens impact param- 
eter. The I class of events with p/uo > 1 is accounted for 
events with shorter time duration and higher magnification 
at maximum. The median values of the two distributions 
are (ti/2)mcd = 1-6 day and (i?max)mod = 20.6 mag. Two 
light curves of I class events are shown in the Figs. [3] and ID 
The first figure (for Mp = 4.74 Mjupitcr) show a more clear 
deviation with respect to the Paczyhski light curve. The sec- 
ond one (for Afp — 0.82 Mjupiter), which is representative 
from a statistical point of view of the whole sample of I class 
events, shows an overall distortion (that in other cases may 
be either symmetric or asymmetric with respect to the max- 
imum) of the hght curve. As far as the II class of events with 
p/uo < 1 is concerned, the dashed lines in the bottom panels 
of the Figs, [l] and [2] show that they have larger time dura- 
tion - (ti/2)mcd ~ 6.4 day - and lower magnification at the 
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Figure 7. Normalized (to unity) distributions of the planet mass 
Mp for the events with detectable planetary deviations (solid line) 
and for the generated events (dashed line). Here we take A^im = 12 
day~^ and D = 8 m. 



Figure 9. Normalized (to unity) time duration ATp niax distri- 
bution of the strongest planet induced perturbation, for I class 
(solid line) and II class (dashed line) events. 
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Figure 8. Upper panel: (normalized to unity) distributions of 
the star-to-planet separation dp (in AU units) for the events with 
detectable planetary deviations (solid line) and for the generated 
events (dashed line). Bottom panel; distributions of d = dp/R^ 
for events as before. 



maximum - (i?max)mcd = 23.1 mag -. Two examples of light 
curves are given in Fig. [5] (for Mp — 0.22 Afjupitcr) and in 
Fig.|6](for Afp — 3.97 Af.jupitcr), with a bump and a multiple- 
peak structure, which is typical of binary microlensing (in 
which the companion mass is large) . T hese features of ca us- 
tic intersections were discussed also bv IPaczviiskil ()l996l ). 

Concerning the reliability of the planetary detections, 
we find that the events of the I class (with p/uo > 1) have 



smaller values of (e)max (for a given Mp value) with respect 
to the II class events. This happens since for the I class 
events the source size p is typically much larger than the 
caustic region, so that averaging the planetary magnifica- 
tion on the source area leads to smaller values of (e)max- 
This does not occurs for the events of the II class (with 
p/uo < 1), for which averaging on the source area is less im- 
portant. This result is reflected in the presence of more clear 
and temporally localized planetary features in the II class 
events. These deviations look similar to that observed in mi- 
crolensing planetary events towards the galactic bulge, for 
which the point-like source approximation is acceptable. We 
also find that (e)max increases with increasing values of Mp, 
a result that is expected since the caustic size is increasing. 

The distributions of the planet mass Mp (for D — 8 
m and the considered lines of sight) are given in the Fig. [7] 
(solid line) for the selected events, (xr > 4, Ngood > 3 and 
{^)rmmax > 0.1). For comparisou, the Afp distribution for 
the whole sample of detectable events (dashed line) is also 
given. From Fig. [7] it follows that larger planetary masses 
lead to higher probability for the detection of planetary fea- 
tures. This result reflects the fact that the detection prob- 
ability is proportional to the caus tic size, which increases 
with the planet-to-star mass ratio llMao fc Paczvriski|[T99ll : 
iBolatto fc Falcdll994l : lGould fc Loeblll992l) . From the same 
figure, it also follows that the planet detection can occur 
with a non negligible probability for Mp > 0.06 A/jupitor 
(Afp > 20 M^), although even Earth mass planets might be 
in principle detectable. However, if we consider telescopes 
with smaller diameter, practically no planet detection oc- 
curs for Afp < 0.06 Afjupitcr and D < 4 m. 

We also recover the well known result that the probabil- 
ity of planet detection is maximized when th e planet-to-star 
separation dp is inside the "len sing zone" (|Gould fc Loebl 
1 19921 : iGriest fc Safizadehl Il998l ). The dp (normalized to 
unity) distribution for selected (solid line) and detectable 
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Figure 10. Histograms for the differences (ti — to) (solid line) and 
(tf — to) (dashed line) for the selected events (xr > 4, Afgood > 3 
and (e)max > 0.1). Initial and final instants for the start and the 
end of the strongest deviation in the light curves are denoted by 
ti and tf, while to is the instant of maximum on the light curve. 



(dashed line) events are shown in the upper panel of Fig. |8l 
The relevance of the lensing zone is clarified in the bottom 
panel of the same figure where the planet separation (in unit 
of the Einstein radius) d — dp /Re is plotted. It turns out 
that about 70% of events with planet detections have d val- 
ues distributed in the lensing zone. We also find an excess 
of I class events at large planetary distances d > 1.6, that is 
related to the interplay between the source size and the size 
of the central caustic. 

The knowledge of the typical time scales for the plan- 
etary perturbations is an important issue to choose an ad- 
equate strategy for the observations, namely, telescope pa- 
rameters and suitable sampling time for optimizing the de- 
tection of the planetary perturbations in the light curves. To 
estimate the time duration of the strongest perturbations we 
introduce a new estimator, Xr n, that is defined as Xr with 
the difference that now the sum runs over the points inside 
the n-th planetary perturbation. We consider a perturba- 
tion to be significative whenever Xr n > 4. The duration 
ATp max is estimated as the time interval with the largest 
value of Xr n- The normalized distribution of ATp max is 
shown in Fig. [O] It results that in pixel-lensing searches to- 
wards M31 typical time duration of the strongest planetary 
perturbations is about 1.5 and 1.2 days for I and II class 
events, respectively. The normalized distributions of the ini- 
tial (ti) and final (tf) instants for the start and the end of 
the strongest planetary deviations (ATp max = tf — ti) given 
in Fig. [To] show that these occur near the maximum mag- 
nification time, as expected since in pixel-lensing the cross- 
ing of the central caustic is more probable. We also find 
that the number of time intervals with significative plane- 
tary deviations on each light curve increases with increas- 
ing values of the ratio p/uo. Indeed, the overall time scale 
ATp tot = ATp n for the significative planetary devi- 



Table 3. As a function of D (first column) we give: the probability 
to detect pixel-lensing events (second column) normalized to the 
events detectable by a 8m telescope, the fraction of I class (third 
column) and II class (fourth column) events, the probability to 
detect planetary features (xr > 4, Ng^od > 3 and (e)max > 0.1) 
for I (fifth column) and II (sixth column) class of events when 
normalized to the events detectable by a telescope with diame- 
ter D and the overall probability (last column). Here we assume 
Mm = 12 day^^ and tcxp = 30 min. 
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ations increases up to 3.4 and 1.6 days for I and II class 
events, respectively. Moreover, our analysis of the distribu- 
tion of ATp tot as a function of telescope size D and sam- 
pling time A^im^ allows us to conclude that a reasonable 
value of the time step for pixel-lensing observations aiming 
to detect planets in M31 is a few hours (Mm — 4 day^^), 
almost irrespectively on D. 

To summarize, the distinctive features of the selected 
events with planetary detections are given by in Table [2] 
(for a telescope with D = 8 m and averaging on the con- 
sidered lines of sight). In particular, we report the median 
values for the distributions of the more relevant quantities 
characterizing the lensing and planetary systems. 

In Table O we present the planet detection probabil- 
ities (by averaging on the selected lines of sight), assum- 
ing Nim ~ 12 day~^, texp ~ 30 min and telescopes of dif- 
ferent diameters For each telescope diameter and class 
of events, the probabilities are evaluated as the ratios be- 
tween the number of the selected events and the number of 
events detectable for the same class and telescope, namely 
Pp = ri(D)/r'{D) and Pf = r'p' {D)/r" (D). The frac- 
tions and /^^ of detectable events for each class is also 
given in Table [3] It results that the probability to detect 
planetary signatures is higher for the events of the I class 
(with p/uo > 1), that however are rare. On the contrary, the 
generated events of the II class are more numerous, but have 
a smaller probability to show detectable planetary features. 
The overall probability (Pp in the last column of the Table 
[3]) is always very small (less than 2 %) and decreases rapidly 
for smaller telescopes. This implies that hundreds of pixel- 
lensing events should be collected to find a few systems with 
planetary features. 

The main result of the present work can be summa- 
rized in Fig. [11] {D = 8 m), where for self-lensing events 
towards M31 with detectable planetary features we present 
the event scatter plot in the (Mp, dp) plane. The thick solid 
lines delimit the region (upper and left) where extrasolar 

Note that, since in pixel-lensing the important parameter is 
the signal-to-noise ratio and it is proportional to D^/texp, to have 
the same probability for planetary feature detection, one can use 
smaller size telescopes as well, by increasing correspondingly the 
exposure time. 



10 G. Ingrosso et al. 



Table 2. Pixel-lensing events with positive planetary detections (xr > 4, A^good > 3 and (e}max > 0.1). Median values of the considered 
distributions. Upper part of the table: I class events [p/uq > 1). Lower part: II class events (p/no < 1). 
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Figure 11. Scatter plot of the planet mass (in unit of Earth 
mass) vs planet distance (in Astronomical Units) . The solid thick 
line delimits the region (upper and left) of planet detection acces- 
sible by radial velocities measurements with a precision up to 1 
m s~^. The observational data were accessed using the extrasolar 
planet on-line catalogue which collects the results of several col- 
laborations (see http://exoplanet.eu/catalog.php and references 
therein). The eight small boxes are the planets detected by the mi- 
crolensing technique. Starting from a sample of 40,000 detectable 
pixel-lensing events {D = 8 m), 630 selected events (indicated 
by black dots) with Xr > 4, Ng^od > 3 and (e)max > 0.1 show 
planetary features and among these 48 events have A/p < 20 . 



Figure 12. The upper panel shows the simulated light curve 
(black dots) of a planetary event with the parameters of the 
be st fit model fo r the PA-99-N2 event (model Wl in Table 1 
of I An et al]|2004l '). In particular, d = 1.84, q = 1.22 X 10"^ (cor- 
responding to a planet mass Mp = 6.34 Mjupiter for a disk lens 
of mass Ml = 0.5 Mq), uq = 3.4 X 10"^, = 132.3 day and 
6 = 24.5 deg. We take the source magnitude A/r = —2.0, and the 
source radius of Rs = 11 Rq (corresponding to p = 1.27 X 10~^ 
and p/uo = 0.37). It is also shown the best fit Paczynski like 
model modified for finite source effects (continuous line), which 
appears almost indistinguishable from the simulated data. The 
bottom panel gives the difference between the two curves. Here 
we use the INT telescope parameters and Ni^^ = 12 day~^. 



planets are detectable by ground based observations, that 
are more sensitive to massive and close-in planets and that 
can be successfully applied only for systems close enough 
to Earth. We remind that current space based observations 
by Keple £3 and COROlH satellites) are expected to de- 
crease the minimum detectable planetary mass limit (up to 
one tenth of the Earth mass) and increase the planetary 
distance (up to tens of AUs). The eight extrasolar planets 
claimed so far to be detected by microlensing since 2003 
in observations towards the Galactic bulge are represented 
by boxes. The locations of points in Fig. [TT] show that the 

http: / /www. nasa.gov/mission_pages/kepler/overview/index.html 
12 http://smsc.cnes.fr/COROT/index.html 



pixel-lensing technique may be used to search for extrasolar 
planets in M31 (including small mass planets), and at the 
moment this is the only method to discover planets in other 
galaxies. As one can see, detectable extrasolar planets have 
planet-to-star separations in the range 0.3 — 25 AU and mass 
in the range 0.1 — 10 Mjupiter (that coincides with the 
assumed lower and upper limits for planetary masses in the 
simulations). However, we note that the detection of plan- 
ets with relative large masses is favourite (see also Fig. [7|. 
We also caution that the planets with Mp < 20 become 
undetectable and disappear from Fig.[TT]if the adopted tele- 
scope has not a good enough photometric stability (about 
0.03 mag, that is the required stability consistent with the 
typical error bars for the detection of small mass planets). 
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Before closing this section we note that an extrasolar 
planet in M31 might have been already detected since an 
anomaly in a p ixel-lensing light curve has been reported 
l|An et alj |2004| ). The authors claim that a binary system 
(lying on the M31 disk) with mass ratio q = 1.22 x 10~^ and 
distance d = 1.84, is a possible explanation of the anomaly 
in the observed light curve. Other parameters are indicated 
in the caption of Fig. 1121 In this figure we give a light curve 
with the best fit parame t ers of the PA-99-N2 event as given 
in Table 1 of lAn et ahl (|2004l ) . It gives a clear deviation 
(Xr = 49, (e)max = 0.6) with respect to the correspond- 
ing Paczyhski shape, at least with our ideal sampling of 
A'ini = 12 day"'^ and observational conditions. In order to 
estimate the secondary object mass, we assume tha t the disk 
star m ass follows the broken power law given bv lAn et al.l 
l|2004l) . Accordingly, one finds a mean mass of ~ 0.5 Mq for 
the lens and therefore a mean value of Mp = 6.34 Mjupitor 
for the planet. This value is at the boundary between the 
planet and brown dwarf region. Our light curve closely re- 
sembles the observed one and the basic characteristics of the 
planetary event fall in the parameter range for the II class 
of events. 



5 CONCLUSIONS 

We consider the possibility to detect planets in M31 by us- 
ing pixel-lensing observations with telescopes of different 
sizes and observational strategies. This is the only way to 
detect planets in other galaxies and acquire information al- 
lowing a comparison of the planetary systems in M31 with 
respect to those in the Milky Way. We carry out MC simu- 
lations and explore the multi-dimensional space of the phys- 
ical parameters of the planetary systems and characterize 
the sample of microlensing events for which the planet de- 
tections are more likely to be observed. We have assumed 
that each lens star in the M31 bulge and disk hosts one 
planet, and used for the planet mass distribution an sim- 
plified law, neglecting any dependence of the planet mass 
on the parent star mass and metallicity. Consideration of 
finite source effects induces a smoothening of the planetary 
deviations with respect to the point-like source approxima- 
tion and, in turn, decreases the chance to detect planets. It 
also implies that in pixel-lensing searches towards M31 only 
few exposures per day could be enough to detect planetary 
features in light curves, at least when using large enough 
telescopes. We find that the pixel-lensing technique favours 
the detection of large mass planets (Mp ~ 2 Mjupitor), even 
if planets with mass less than 20 can be detected (with 
small probability, however) by using large telescopes with 
a sufficient photometric stability. Microlensing is intrinsi- 
cally a "no repetition" phenomenon and variable stars may 
mimic microlensing events and contaminate the sample of 
events attributed to microlensing. Therefore, real observa- 
tions should be done at least in two bands, to check for 
achromaticity and be confident that the contamination by 
variables can be sorted out. However, a minor chromaticity is 
expected since the source limb darkening profile depends on 
the considered band and on the spectral type of th e source 
star (see, e. g., iBogdanov fc Cherepashchuklll995bl : iPeichal 
l2009h . 

Finally, we remark that although we have neglected the 



contribution to microlensing events of MACHOs in both 
galactic halos (in this respect the estimated planet rate 
should be considered as a lower bound), pixel-lensing ob- 
servations towards M31 could be very useful in establishing 
whether planets may form around MACHOs as well. 
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